Quantum state reduction for universal measurement based computation 
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Measurement based quantum computation (MBQC), which requires only single particle measurements on a 
universal resource state to achieve the full power of quantum computing, has been recognized as one of the 
most promising models for the physical realization of quantum computers. Despite considerable progress in 
the last decade, it remains a great challenge to search for new universal resource states with naturally occurring 
Hamiltonians, and to better understand the entanglement structure of these kinds of states. Here we show that 
most of the resource states currently known can be reduced to the cluster state, the first known universal resource 
state, via adaptive local measurements at a constant cost. This new quantum state reduction scheme provides 
simpler proofs of universality of resource states and opens up plenty of space to the search of new resource 
states, including an example based on the one-parameter deformation of the AKLT state studied in [Commun. 
Math. Phys. 144, 443 (1992)] by M. Fannes et al. about twenty years ago. 



PACS numbers: 03.67.Lx, 03.67.Pp 

Measurement based quantum computation (MBQC) [ 1 ], an 
interesting computation model that incorporates peculiar as- 
pects in quantum mechanics like entanglement and measure- 
ment, achieves the full power of quantum computing by adap- 
tive local measurements on a resource state. The first MBQC 
scheme, also known as the one-way quantum computation, 
employs the now well-known cluster state [2], The highly 
entangled feature of the cluster state indicates that high en- 
tanglement is a key requirement for universality in MBQC. 
Although this is true in some sense [3, 4], it is also clear now 
that too much entanglement could also undermine universal- 
ity [5, 6]. In other words, the entanglement should be man- 
ageable in a structured way. 

As shown in the recent breakthrough made by D. Gross et 
al. in Refs. [7, 8], the matrix product state (MPS) formal- 
ism [9, 10] or, in higher spacial dimensions, the computa- 
tional tensor networks (CTN) [8, 11, 12] provides such an 
infrastructure for manipulating the entanglement and brings 
a new scheme of MBQC called the correlation space quantum 
computation. In this framework, a lot of new resource states 
beyond the cluster state are proposed. Most of the new re- 
source states have different properties from the cluster state 
concerning, for example, local entropy, the two-point corre- 
lation function and the locality of the Hamiltonians of which 
they are unique ground states. Especially, some of the new 
resources are unique ground states of more practical Hamilto- 
nians [13, 14], thereby overcoming the major flaw of the clus- 
ter state of not being a unique ground state of any two-body 
nearest-neighbour gapped Hamiltonian [15]. 

Here, we introduce the concept of quantum state reduc- 
tion for MBQC and the motivation is twofold. First of all, 



state reduction serves as a tool for revealing the entanglement 
structure of universal resource states [16, 17]. Similar to the 
common technique to study entanglement by considering lo- 
cal transformations [18-20], quantum state reduction is also 
a type of local transformation tailored to meet the nature of 
MBQC. We find out that almost all known resource states can 
be locally transformed to a cluster state via quantum state re- 
duction, indicating that these resource states possess similar 
entanglement structure as the cluster state. Compared to the 
attempt made in Ref. [2 1 ], where it was shown that almost all 
these resource states are "universal state preparators", quan- 
tum state reduction is a more direct approach and most notably 
more respectful for the geometry of the resource. 

Secondly, although the application of the MPS/CTN for- 
malism in the theory of MBQC is elegant and fruitful, the 
routine for analyzing the universality of a resource state re- 
mains a complicated procedure, including the initialization, 
embedding of universal rotations, the readout, and compen- 
sation for the randomness. The quantum state reduction ap- 
proach largely simplifies the analysis. For example, the uni- 
versality of AKLT state [22, 23] is now cleanly summarized 
in Fig. 2. The simplicity also enables us to find new resource 
states, giving the universality of two deformations of AKLT 
state almost for free. 

To be more precise, our state reduction is a transformation 
from one resource state to some other universal target state 
(usually the cluster state) using local measurement and adap- 
tive classical control. This transformation is named reduction 
as it resembles the reduction in complexity theory — as long as 
|\&) is reducible to |$), it is in principle no harder to construct 
MBQC schemes for |^) than for |$). It is important to note 
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that, although the reduction is a random procedure, it always 
succeeds in obtaining a target state. Possibly, the resulting 
state may consist much smaller number of particles than the 
state before reduction. However, the cost or efficiency of the 
reduction, measured by the diminution in the number of par- 
ticles, is always expected to be a constant. 

Matrix product states. — Following the notion in Refs. [7, 
8], a matrix product state | *!>„) of n particles has the following 
form 

|tf n }= {R\A[x n ]---A[x 1 ]\L)\x l ---x n ). (1) 

xi , ■ ■ ■ jXn—O 

The physical dimension of each site is d, while the bond di- 
mension of the state defined by the size of the matrices is 8. 
In general, one may consider an MPS where the defining ma- 
trices are site-dependent. The defining matrices are usually 
far more important than the boundary conditions (R\, \L) and 
sometimes we will use only the d-tuple (A[0], . . . , A[d — 1]) 
to specify an MPS. We will ignore the effect of a local change 
of basis on the physical space. For example, a state may be 
said to have some MPS representation even though this is only 
true up to some local unitary operations. Also for simplicity, 
A[x] 's will be given up to some normalization constant, which 
can be figured out from J2 X Alx^Alx] = I shown in Ref. [ ]. 

A lot of states of particular interest in quantum information 
are indeed matrix product states of small bond dimension. The 
cluster state of one spacial dimension (i.e. a chain, or a 1-D 
cluster state), for example, is an MPS with defining matrices 
(H, HZ), where H, X, Y, Z are used to denote the Hadamard 
and Pauli matrices respectively. The GHZ state can be rep- 
resented by (I, Z) up to local rotations. Another intriguing 
example is the AKLT state [22, 23] first studied in condensed 
matter theory. It will be shown in the next section that it has a 
simple MPS representation (I, X, Z). 

The correlation space quantum computation employs the 
structure of MPS as in Eq. (1). It starts with the initial state 
\L) of the so called correlation space, measures the physi- 
cal spaces sequentially and thereby processes the correlation 
space, and finally reads out the information stored in the corre- 
lation space. Several new universal resources for MBQC were 
introduced in Refs. [7, 8] including a modified AKLT state 
with defining matrices (H, X, Y). Later, the original AKLT 
state is also shown to be universal for MBQC [ 1 3]. Recently, 
the concept of quantum wires is defined and fully character- 
ized in Ref. [24], which essentially gives the explicit condition 
for an MPS with d = 5 = 2 to be a universal resource. There 
are two normal forms of the matrices for quantum wires. One 
is the "byproduct normal form" 

A[0] = W/V2, A[l] = WS(<f>)/V2, (2) 

and the other is the "biased normal form" [25] 

A[0] = siwyW', A[l] = cos 7 W'Z, (3) 

where W, W are rotations along axes in the X-Z plane of the 
Bloch sphere and S((j)) = exp(— i(f>Z/2). 



There is a higher spacial dimensional generalization of 
MPS, known by the names of the computational tensor net- 
works [8], or the projected entangled pair state (PEPS) [11, 
12]. Due to limitations of space, we refer the readers to the 
references above for details. 

The tabular form. — For the convenience of later discus- 
sions, we introduce a tabular form of MPS. In the tabular form, 
one writes the defining matrices of a block of sites explicitly 
in a table, where each column consists the d matrices of a cor- 
responding site. The physical indexes determine a selection of 
one matrix from each column, whose product gives the correct 
amplitude together with the boundary conditions (R\ and \L). 
From the definition in Eq. ( 1) and the properties of MPS [ ], 
we have 1) For any two neighboring columns, multiplication 
of M to the right of all matrices in the left column and M _1 
to the left of all matrices in the right column simultaneously 
does not change the state; 2) A unitary transformation in the 
physical space corresponds to linear combinations of entries 
in the column with coefficients of the unitary; 3) Measure- 
ment in the computational basis corresponds to the deletion 
of column entries not consistent with the measured outcome; 
and 4) Columns of single entry can be removed by absorbing 
them to a neighboring column. We will use binary relations 
=, ~ and -< to represent equality, local unitary equivalence 
and quantum state reduction respectively. 

As an example, Table 1 of Fig. 1 consists of a block of two 
sites of the AKLT states which will be discussed in the next 
section and it equals Table 2 by Property 1) of the tabular 
form. 

Reduction of the AKLT state.— The AKLT state [22, 23], 
named after Affleck, Kennedy, Lieb, and Tasaki, has become 
one of the prototypical states of spin systems. It also gives an 
excellent example for quantum state reduction. 

As the origin of matrix product states, the AKLT state bears 
a simple MPS representation with 

A[0] = Z,A[1] = x/2|0)(lU[2] = V2\l)(0\. (4) 

Up to a local unitary operation, the matrices of the AKLT state 
can also be chosen as (X, Y, Z). In fact, any three different 
matrices of the identity and the Pauli matrices will work. For 
example, Fig. 1 presents the proof that (I, X, Z) also stands 
for the AKLT state. In this figure, Table 2 is obtained by 
adding the Y's with blue color, and hence represents the same 
state as Table 1 by Property 1) of tabular form; Table 2 and 
3 describe two states that are equivalent under local unitary 
operations by Property 2). 
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FIG. 1. AKLT as (I,X, Z) 



We now show the reduction from the AKLT state to the 1-D 
cluster state. It is convenient to start with the (I, X, Z) form. 
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Two different measurements Ai and A/2 will be used alterna- 
tively, where M\ measures {|0), |1)} versus |2), and A/2 mea- 
sures {|0), 1 2)} versus |1). That is, each measurement con- 
sists a two dimensional and a one dimensional projectors. The 
measurements are called success (failure) if the outcome cor- 
responds to the two (one) dimensional subspaces. We measure 
the two measurements sequentially along the AKLT chain and 
switch the measurement we use only when the previous one 
succeeds. This simple procedure is called the alternating mea- 
surement scheme. 

Table 1 in Fig. 2 denotes a possible result after the alternat- 
ing measurements on the AKLT state. More specifically, one 
first measures A/i and succeeds. Next, the measurement A2 is 
used. It results in the single dimensional space |1) once, and 
succeeds subsequently, and so on. After renaming the phys- 
ical indexes and absorbing the X and Z in red color to their 
previous columns, we have Table 2 in Fig. 2 by Property 4) 
and 2). This is actually already a 1-D cluster state by the sec- 
ond line of reasoning in Fig. 2. 
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FIG. 2. AKLT reduced to 1-D cluster 



Families of universal states of AKLT type. — The simplicity 
of the above analysis enables us to generalize the same ap- 
proach to a larger family of AKLT type of states. Notice that 
the key property that validates the first line of Fig. 2 is simply 
X 2 = Z 2 = I and the key to the second line is that H 2 = I 
and XH = HZ. We now choose two unitary matrices A 
and B such that A 2 = B 2 = I, where A, B correspond to 
7r-rotations along n a and nj, on the Bloch sphere respectively. 
Let C oc A + B be the it -rotation along n a + n&. We will 
have C 2 = I and AC = CB. Therefore, we can prove the 
following reduction similarly 

(I,A,B) <{C,CB). (5) 

Employing the gauge freedom of the representation of MPS, 
one can always choose B to be Z and C to be sin 9X+cos 9Z, 
making (C, CB) a quantum wire in the normal form of 
Eq. (2). Note that the 1-D AKLT state is a special case where 
6 = 7r/4 and that the error group (C,B) is isomorphic to 
the dihedral group for infinitely many 0's. With techniques of 
Ref. [9], one can check that the new AKLT type resource is 
always unique ground state of a nearest-neighbor, frustration- 
free Hamiltonian. 

In Ref. [10], Fannes, Nachtergaele and Werner consid- 
ered another one-parameter deformation of the AKLT model 



whose ground state is an MPS with 

A[0] = BinOZ,A[i\ = cos0|O)(l|,A[2] = cos0|l)(O|. 

We will show the universality of these states also by reduction. 
First, the defining matrices can be chosen as 

(sin 0Z, cos 8X/V2, cos 6Y/V2) , 

up to local unitary transformation. Let 8 be the angle that sat- 
isfies tan0 = \/2tan0. The defining matrices can be simpli- 
fied to (sin 91, cos OX, cosdZ), in the same way as in Fig. 1. 
Using a similar alternating measurement scheme in Fig. 2, this 
is further reducible to (sin 8H, cos OHZ), a universal state in 
the biased normal form of Eq. (3). One caveat is that, in this 
case, we cannot simply absorb the X and Z in red color in 
Fig. 2 to neighboring sites because of the bias. But one can 
always measure the computational basis in several neighbor- 
ing sites and cancel their effects by a random walk on the Pauli 
group. 

Reduction of quantum wires to cluster states. — We now dis- 
cuss the reduction of universal quantum wires to the clus- 
ter state. The special case of (W, WZ) is much easier to 
deal with. To transform it into (H, HZ), one can simply 
implement HW> in the correlation space using the sites be- 
forehand. In the general case of (W, WS ((/>)), however, one 
cannot succeed using projective measurement only — the local 
entropy determined by (j> [24] can never be increased. Yet, 
if the more general quantum measurement is employed, this 
is again possible. It will be easier to work with the biased 
normal form in this case. Suppose we want to transform 
{sm^W iCOS^W Z) to (H,HZ). Assume that 7 g (0,7r/4] 
with out loss of generality and apply on the site a general mea- 
surement with operators 

Mo = |0)(0|+tan 7 |l)(l|,M 1 = yjl - tan 2 7 |1)(1|, 

known as the filtering operation. When the outcome happen 
to be 0, we have changed the matrices to (W, WZ) and can 
proceed as in the easy case; otherwise, we need to undo the 
action of W'Z on the correlation space and start all over again. 
Note that it's also possible to reduce a universal quantum wire 
to another quantum wire that is different from the cluster state 
similarly. 

Higher spacial dimensional cases. — This section investi- 
gates the idea of quantum state reduction in the case of higher 
dimensional resource states, which are necessary for the full 
power of universal quantum computing. 

The triCluster state, an interesting variant of the cluster 
state, is proposed in Ref. [14] as a universal resource state 
of local dimension 6 and is the unique ground state of a two 
body, frustration-free, gapped Hamiltonian. It's not difficult 
to see that there is a reduction to cluster state on exactly the 
same lattice of the triCluster state and we leave the details to 
Appendix A. 

For most of the known 2-D resource state, a general cou- 
pling scheme has been used to make 2-D resource states out 
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of 1-D chains as in constructing the 2-D AKLT resource, ana- 
lyzing 2-D weighted graph state [7, 8], and weaving quantum 
wires into quantum webs [24]. In this scheme, one can always 
(a) isolate several usually horizontal, 1-D, universal chains 
from the 2-D state and (b) couple the correlation space of 
two neighboring 1-D chains whenever necessary. Resources 
of this type can be transformed to 2-D cluster state. To see 
this, one first isolate 1-D chain states from it; use the reduc- 
tions we already have for the 1-D case to obtain 1-D cluster 
states; and then employ an appropriate coupling to link the 
1-D cluster states into a two dimensional cluster state. The 
first two steps are obvious, while the third step is possible as 
shown in Fig. 3. In this figure, A[x]'s are defining matrices of 
cluster state chosen to be A [0] = |+)(0|,v4[l] = |-)(1|, and 
B, C are tensors in the notion of Refs. [7, 8] 



B[l} = 
C[l] = 



|->r(l|l®|->* 
|-)r(l|l® |l)u- 



B[0] = \+) r (0\l 8> \+)d, 

C[Q] = \+)r(Q\l 8> |0> ttJ 



The Hadamard gates and the CZ gate |0) (0| (g) 1+ |1)(1| (8 Z, 
are implemented on the correlation spaces. The right hand 
side of Fig. 3 represents two nodes of degree 3 of a cluster 
state. 



1-D state. And in the 2-D case, the preparator method will 
need polynomial cost to obtain a 2-D cluster state while state 
reduction remains of constant cost. 

The tabular form we propose is well hinged to the structure 
and properties of MPS. It simplifies the analysis by hiding the 
unwanted details and provides an intuitive way of manipulat- 
ing the matrices. We have mainly investigated reductions of 
MPS and PEPS resources, but the idea seems to be able to gen- 
eralize to potentially new MBQC scheme not known yet. It's 
also reasonable to believe that investigations of the reduction 
method will improve our understanding of both the MBQC 
itself and the structure of universal resource states. 
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ogy, Sydney, and the NSF of China (Grant Nos. 60736011 
and 60702080). ZJ acknowledges support from NSF of China 
(Grant Nos. 60736011 and 60721061). BZ is supported by 
NSERC and Quantum Works. Research at Perimeter Institute 
is supported by the Government of Canada through Industry 
Canada and by the Province of Ontario through the Ministry 
of Research & Innovation. 
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FIG. 3. Reduction of 2D resources of the coupling scheme 

In the weighted graph state, for example, the isolated wire 
has the form HZ X S Z where z depends on the outcomes of 
neighboring sites [7, 8]. We can measure all sites with odd z 
in the 0-1 basis and the resulting state is a 1-D cluster state up 
to random Clifford byproduct operations. The CZ gate can be 
applied in the same way as in Ref. [8]. Other coupling based 
schemes mentioned above can be analyzed in a similar way. 

Discussions. — The method of reduction for proving univer- 
sality of MBQC resource state is applicable to other examples 
that have not been covered in the previous sections. These in- 
clude, for example, the modified AKLT state (H, X, Y) pro- 
posed in Ref. [7] and the second toric code state example in 
Ref. [8]. At the current stage, however, we do not know how 
one can reduce the first toric code example to the 2-D cluster 
state. 

It is worth comparing the idea of reduction and that of the 
"universal state preparator" result of Ref. [21]. A reduction 
to cluster state would imply the universal preparator property 
of the resource. On the other hand, although any universal 
preparator could in principle be transformed to a cluster state, 
the transformation does not respect the underlying lattice of 
the resource and may be less efficient than the quantum state 
reduction we are considering. For example, in the 1-D case, 
the method in Ref. [2 1 ] isolates a single particle of the phys- 
ical space while reduction of 1-D resource results in again a 
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Appendix A: Reduction of triCluster state 

The triCluster state considered in Ref. [14] is a universal 
resource state of local dimension 6 and is the unique ground 
state of a two body, frustration-free, gapped Hamiltonian. We 
give a detailed reduction from the triCluster to cluster state in 
this appendix. 

The most concise way of describing the triCluster state is 
to employ the PEPS picture [12]. As in Fig. 4, each bond is 
the state \H) oc |+)|0) + |— and each dashed circle is the 
projection P = P + P x + P 2 where |±) cx |0) ± |1) and 

P = |0)(000| + |1)(111|, 
P = |2)(001| + |3)(110|, 
P 2 = |4)(010| + |5)(101|. 




FIG. 4. Reduction of triCluster state to cluster state 

One can transform the triCluster state into a cluster state on 
the same lattice, by simply measuring each site with Qj = 
PjPj for j = 0, 1, 2. Although the outcomes will be random, 
one can assume that we have always measured 0, thereby only 
Pq are used for the projection on each site, except that there 
will be some random X errors happening on the bonds be- 
fore the application of projections. As the action of X on one 
end of the bond \H) is equivalent to a Z on the other end, we 
can propagate all X's to neighboring sites as Z's on the bond, 
which is the same as Z's on the physical space. In conclu- 
sion, we have obtained the cluster state up to some Z errors 
determined by the random outcomes. 



Appendix B: On the synchronization problem in the coupling 
method 

We have discussed the reduction for a large class of 2-D re- 
sources made out of 1-D state by the coupling method. There 
is, however, a tricky point that we have overlooked. This ap- 
pendix aims to convey the general idea that the problem is 
solvable. 

Imagine that we first isolate two neighboring chains, re- 
duce them to 1-D cluster separately and then try to connect 
them into a 2-D cluster. Because of the randomness of the 
measurements, it may happen that the two chains are not both 
ready for the interaction if they are not aligned to the same 
column. That is, one chain contains more unmeasured par- 
ticles than the other. To solve this synchronization problem, 
one may measure along these two chains in the basis that will 
induce the error-group random walk. Let us consider a simple 
example where the error group is (H, Z) and the interaction 
implemented is the CZ gate. The random walk is the Markov 
chain on the graph depicted in Fig. 5. The Markov chain is 
of period 2, meaning that it will only return to the starting 
state after even number of steps. Now the measurements on 
the two chains induce two independent random walk. If the 
difference of the number of unmeasured particles is even, the 
two random walks will be in the state / simultaneously after 
a finite number of steps. Otherwise, they will never be I si- 
multaneously because of the periodicity of the Markov chain. 
However, one can wait a finite steps for the configuration of / 
and Z. The Z operation can commute with the CZ interaction 
and amounts to a byproduct for future steps of the reduction. 




FIG. 5. Random walk on the error group (H, Z) 

A similar argument will work for the case where Z is a 
generator of the error group and the interaction implemented 
is CZ. As Z is a generator, the period of the random walk is 
either 1 or 2, and the synchronization problem can be dealt 
with similarly. 



